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ABSTRACT 

Our analysis of high-resolution Goddard High-Resolution Spectrograph (GHRS) spectra of late-type 
stars shows that the Si iv and C iv lines formed near 10 5 K can be decomposed into the sum of two 
Gaussians, a broad component and a narrow component. We find that the flux contribution of the 
broad components is correlated with both the C iv and X-ray surface fluxes. For main-sequence stars, 
the widths of the narrow components suggest subsonic nonthermal velocities, and there appears to be a 
tight correlation between these nonthermal velocities and stellar surface gravity (£ NC cc g °- 68±0 - 07 ). For 
evolved stars with lower surface gravities, the nonthermal velocities suggested by the narrow components 
are at or just above the sound speed. Nonthermal velocities computed from the widths of the broad 
components are always highly supersonic. We propose that the broad components are diagnostics for 
microflare heating. Turbulent dissipation and Alfven waves are both viable candidates for the narrow 
component heating mechanism. 

A solar analog for the broad components might be the “explosive events” detected by the High- 
Resolution Telescope and Spectrograph (HRTS) experiment. The broad component we observe for the 
Si iv A1394 line of a Cen A, a star that is nearly identical to the Sun, has a FWHM of .109 ± 10 km s 1 
and is blueshifted by 9 ± 3 km s _1 relative to the narrow component. Both of these properties are con- 
sistent with the properties of the solar explosive events. However, the a Cen A broad component 
accounts for 25% ± 4% of the total Si iv line flux, while solar explosive events are currently thought to 
account for no more than 5% of the Sun’s total transition region emission. This discrepancy must be 
resolved before the connection between broad components and explosive events can be positively estab- 
lished. 

In addition to our analysis of the Si iv and C iv lines of many stars, we also provide a more thorough 
analysis of all of the available GHRS data for a Cen A (G2 Y) and a Cen B (K1 V). We find that the 
transition region lines of both stars have redshifts almost identical to those observed on the Sun. 
showing an increase with line formation temperature up to about log T = 5.2 and then a rapid decrease. 
Using the O rv] lines as density diagnostics, we compute electron densities of log n e = 9.65 ± 0.20 and 
log n e = 9.50 ± 0.30 for a Cen A and a Cen B, respectively. 

Subject headings: line: profiles — stars: activity — stars: chromospheres — stars: coronae — 
stars : late-type — ultraviolet : stars 


1. INTRODUCTION 

Edlen’s (1941) identification of highly ionized species in 
the solar corona and the subsequent detection of solar and 
stellar X-rays and ultraviolet emission lines demonstrated 
that the outer atmospheres of late-type stars are heated to 
temperatures greatly in excess of the photospheric surface 
temperatures. These heated regions are called chromo- 
spheres (for the T < 2 x 10 4 K plasma), transition regions 
(for the 2 x 10 4 K < T < 10 6 K plasma), and coronae (for 
the T > 10 6 K plasma). The existence of the hot gas 
requires nonradiative, or “mechanical,” heating processes. 
Despite five decades of effort, there has been only limited 
success in deducing the nature of these heating processes. 
Proposals include the dissipation of acoustic waves via 

1 Based on observations with the NASA/ESA Hubble Space Telescope, 

obtained at the Space Telescope Science Institute, which is operated by the 

Association of Universities for Research in Astronomy Inc., under NASA 

Contract NAS5-26555. 


shocks (proposed as early as the 1940s), turbulent or shock 
wave dissipation of various modes of MHD waves, Joule 
heating by electric currents, quasi-steady magnetic field 
reconnection, and rapid field reconnection either as low- 
energy transient events (microflaring) or high-energy events 
(flaring). On the Sun, the transition region and coronal 
emission arises predominantly in regions dominated by 
magnetic fields, which suggests the existence of a magnetic 
heating mechanism. For a comprehensive review of possible 
heating mechanisms, see Narain & Ulmschneider (1990, 
1996) and the papers from the Heidelberg conference 
“Mechanisms of Chromospheric and Coronal Heating” 
(Ulmschneider, Priest, & Rosner 1991). 

The main difficulty in identifying which mechanism(s) 
heats the chromosphere or corona of a specific late-type star 
or a class of such stars is that one observes only the gross 
effects of the heating process, with no spatial resolution and 
(usually) poor time resolution. Stellar X-ray luminosities 
and the broadband spectral energy distributions observed 
with the Einstein, ROSAT, and ASCA satellites, and line 
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Star 

Spectral 

Type 

Distance 

(pc) 

Mass 

(M e ) 

Radius 

(*e) 

T e{( 

(K) 

v sin i 
(km s -1 ) 

”rad 

(km s _1 ) 

References 

AU Mic 

M0 Ve 

9.4 

0.58 

0.57 

3770 

<5.9 

-2.1 

1, 2, 3, 4, 5 

Capella 

G1 III (G8 III) 

13.3 

2.56 (2.69) 

9.2 (12.2) 

5700 (4940) 

36 (3) 

Variable 

6, 7 

Procyon 

F5 IV-V 

3.5 

1.75 

2.1 

6500 

2.8 

-3.5 

8, 9, 10 

HR 1099 

K1 IV (G5 IV) 

36 

1.4 (1.1) 

3.9 (1.3) 

4490 (5733) 

38 (13) 

Variable 

11, 12 

a Cen A 

G2 V 

1.3 

1.085 

1.24 

5770 

<2.2 

-23.4 

1, 13, 14, 15, 16, 17 

a Cen B 

K1 V 

1.3 

0.90 

0.84 

5300 

<1.0 

-20.3 

1, 14, 15, 16, 17, 18 

31 Com 

GO III 

-90 

-2.6 

9.1 

5700 

57 

-1 

19, 20, 21, 22, 23 

P Cet 

K0 III 

16.4 

2.5 

8.9 

4900 

3 

+ 13.1 

1, 24, 25, 26, 27, 28 

j? Dra 

G2 Ib-II 

-200 

7.2 

70 

5275 

13 

-20.3 

29, 30, 31, 32 

P Gem 

K0 III 

10.0 

1.7 

9.1 

4865 

2.5 

+ 3.3 

1, 24, 25, 33, 34 

AB Dor 

K0-K2 IV-V 

15.3 

-0.9 

-1.3 

-5460 

90 

-30 

35, 36, 37, 38, 39 

Sun 

G2 V 


1.0 

1.0 

5800 

2.0 




References. — (1) Gliese & Jahreiss 1991; (2) Pettersen 1980; (3) Bopp & Fekel 1977; (4) Henry & McCarthy 1993; (5) Linsky & Wood 1994; (6) 
Hummel et al. 1994; (7) Strassmeier & Fekel 1990; (8) Irwin et al. 1992; (9) Gray 1982a; (10) Steffen 1985; (11) Fekel 1983; (12) Vogt & Penrod 1983; (13) 
Hallam, Altner, & Endal 1991 ; (14) Heintz 1982; (15) Lydon, Fox, & Sofia 1993; (16) Demarque, Guenther, & van Altena 1986; (17) Noels et al. 1991 ; (18) 
Char et al. 1993 ; (19) Ayres et al. 1995; (20) Strassmeier et al. 1990; (21) Hirshfeld, Sinnott, & Ochsenbein 1991 ; (22) Oranje, Zwaan, & Middlekoop 1982; 
(23) Flower 1977; (24) Gray 1982b; (25) Astronomical Almanac 1996; (26) Eriksson, Linsky, & Simon 1983; (27) Iben 1967 ; (28) Bell & Gustafsson 1978; 
(29) Brown et al. 1984; (30) Hoffleit & Jaschek 1982; (31) Beavers & Eitter 1986; (32) Luck 1982; (33) Blackwell et al. 1990; (34) Drake & Smith 1991; (35) 
Rucinski 1985;(36)Guirado etal. 1996; (37) Vilhuetal. 1993 ; (38) Vilhuet al. 1996;(39) Hempelmann et al. 1995. 


fluxes observed with EUVE, provide information on the 
corona’s response to heating, but no clear signatures of 
the nature of the heating mechanism(s) itself. Similarly, the 
fluxes of ultraviolet emission lines observed with IUE and 
Hubble Space Telescope ( HST ) provide information on the 
emission measure distribution in the chromosphere and 
transition region, which again is simply a response to the 
heating. What is needed are observational discriminators 
that can identify the different heating mechanisms. 

The intensity of the observed emission and its correlation 
with stellar parameters provides an initial crude discrimi- 
nator between different types of heating processes. For 
example, Schrijver (1987) and others have noted that the 
radiative losses per unit area of a star (the stellar surface 
flux) in important emission lines formed in chromospheres 
and transition regions have minimum values that depend 
upon B — V color (or ^eff) and upon gravity. This so called 
“basal flux,” observed from old slowly rotating stars and 
also seen in the centers of solar supergranule cells away 
from strong magnetic fields, is generally thought to be the 
result of heating by shocks from simple acoustic waves gen- 
erated in the stellar convective zone (see, e.g., Ulmschneider 
1990). Such heating must always be present and is the 
logical explanation for the weak background emission level 
seen in the low-activity stars. 

Enhanced X-ray and UY-line surface fluxes characterize 
“active stars,” a generic term that includes pre-main- 
sequence stars, young main-sequence stars, dMe and other 
flare stars, and rapidly rotating binaries (e.g., RS CVn, 
Algol, and W UMa systems). The active stars have X-ray 
and UV-line surface fluxes up to 10 4 times larger than the 
basal flux stars. For some diagnostics like X-rays, Mg n, 
and C iv, there appears to be an empirical upper limit to the 
surface flux that is now called “ saturated heating ” (Vilhu 
1987). The excess heating above the basal rate is empirically 
correlated with stellar rotation and is generally thought to 
be magnetic in character, but until now it has been unclear 
whether the heating is due to the dissipation by turbulence 
or shock waves in a magnetic field (e.g., MHD waves; Stein 
1981) or to the rapid annihilation of magnetic fields to 
produce accelerated particles and subsequent heating (see, 


e.g., Parker 1988), a process often called “ microflaring.” In 
this paper, we describe a simple diagnostic signature that 
might allow the relative importance of these two types of 
magnetic heating processes to be estimated empirically. 

Linsky & Wood (1994) discovered that HST/GHRS pro- 
files of the Si iv and C iv emission lines of the dMOe flare 
star AU Mic, which are formed at a temperature of about 
10 5 K in the transition region, contain broad wings extend- 
ing out to at least ±200 km s -1 . They found that the AU 
Mic profiles could be fitted with two Gaussians, a narrow 
component and a broad component. In Wood’s (1996) rea- 
nalysis of these data, the widths of the two components 
were found to be FWHM = 31 ± 3 km s -1 and 
FWHM = 162 + 11 km s~ 1 , which suggests most probable 
nonthermal velocities of £ = 15.4 + 1.5 km s~ 1 and 
£ = 97 + 9 km s" 1 , respectively. The nonthermal velocity 
of the narrow component is somewhat less than nonthermal 
velocities observed in the solar transition region (£ 0 « 26 
km s -1 ; Dere & Mason 1993), but it is the nonthermal 
velocity of the broad component that is strikingly different. 
Linsky & Wood (1994) proposed that microflaring in the 
transition region of AU Mic is responsible for the broad 
components, which implies that microflaring is an impor- 
tant source of atmospheric heating for AU Mic. 

Linsky et al. (1995) and Wood (1996) also detected broad 
components in the transition region lines of Capella and 
showed that the broad wings are not present in chromo- 
spheric lines such as the Mg ii h and k lines and the O i 
A 1300 triplet. For the G1 III star, which provides about 
91% of the flux in the transition region lines, they found 
that the broad components are responsible for about 41% 
of the total Si iv and C iv flux. Broad components were also 
detected in the spectra of the cooler star (K1 IV) of the RS 
CVn system HR 1099 (Wood et al. 1996). The chromo- 
spheric Mg ii h and k lines of HR 1099 also have broad 
wings, which suggests that for HR 1099, the broad com- 
ponent phenomenon might extend into the chromosphere, 
although the very high opacity of the Mg ii lines makes this 
interpretation uncertain. The C rv broad components of the 
K1 IV star, which produces essentially all of the line flux of 
HR 1099, contribute about 66% of the total line flux. In 
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contrast, Wood et al. (1996) found no broad components 
( <20% of the total line flux) in the C iv and Si iv profiles of 
the relatively inactive star Procyon (F5 IV-V); although the 
N v line, which is formed at higher temperatures, does 
appear to have a broad component that accounts for about 
42% of the N v flux. 

In this paper, we use GHRS observations of seven addi- 
tional stars to investigate stellar transition region line pro- 
files further, specifically to explore the nature of the broad 
wing emission seen in the C iv and Si iv lines of AU Mic, 
Capella, and HR 1099. The full sample of stars is sum- 
marized in Table 1. Two sets of stellar parameters are listed 
for Capella and HR 1099, as both are close, unresolved 
binary systems. The parameters given in parentheses are for 
the star that is only a minor contributor to the UV emission 
lines. Vilhu et al. (1996) used two-Gaussian fits to model the 
C iv lines of the Pleiades age star AB Dor (K0-K2 IV-V), 
similar to our analyses of the C iv lines of AU Mic, Capella, 
and HR 1099, so those results also are considered here. We 
present a complete analysis of the available GHRS observ- 
ations of a Centauri A (G2 V) and a Centauri B (K1 V). We 
also discuss GHRS observations of j8 Ceti (K0 III), 31 
Comae (GO III), /? Geminorum (K0 III), and /? Draconis (G2 
lb— II). However, because the GHRS data sets for these four 
stars are considerably larger than the very limited a Cen 
A + B data set, we confine our attention to only the C iv and 
Si iv lines. More thorough analyses of the Cet, 31 Com, 
Gem, and j8 Dra data will be presented in future papers. 

2. GHRS OBSERVATIONS 

The GHRS instrument and its capabilities prior to the 
installation of the COSTAR corrective optics in 1993 
December are described in detail by Brandt et al. (1994) and 
Heap et al. (1995). COSTAR significantly improved the 
spectral resolution of observations made through the large 
science aperture (LSA). Table 2 lists the GHRS spectra 
analyzed here, most of which are post-COSTAR, LSA 
observations. 

The fundamental data are moderate-resolution observ- 
ations of the C iv 21550 and Si iv 21400 spectral regions, 
except for a Cen A and a Cen B for which we lack the C rv 
spectra. For Dra, we have C iv spectra obtained in three 
different years. Echelle spectra of a Cen A and a Cen B 
containing the Lyman-a, O v] 21218, and Mg n h and k 


lines are listed also in Table 2. These observations already 
have been presented by Linsky & Wood (1996) in an 
analysis of interstellar absorption features. We refer the 
reader to that paper for a description of the data reduction 
performed on those observations. The procedures used for 
the moderate-resolution data are discussed briefly below; 
they are described in more detail in Wood et al. (1996), 
Wood (1996), Linsky et al. (1993), and references therein. 

The Pt-Ne calibration lamp was observed before each of 
the observations listed in Table 2 (except for the a Cen A 
G140M observation — see § 3.1), and these spectra were used 
to calibrate the wavelengths of the science spectra. All of the 
observations were subdivided into two or more readouts 
(see Table 2). Breaking up the observations in this fashion 
allows us to correct for line shifts induced during the course 
of the observations by changes in Earth’s magnetic field 
and/or instrumental thermal effects. The wavelengths of the 
individual readouts are corrected by cross-correlating them 
against the first readout, which is taken nearest in time to 
the wavelength calibration image. We also searched the 
individual readouts for flux variations that might indicate 
flaring. No such variations were detected, however. 

The /? Cet, 31 Com, and /? Gem spectra were obtained 
using the FP-SPLIT mode, in which the spectrum is shifted 
on the diode array between readouts several times during 
the observation. This technique allows one to compensate 
for sensitivity variations in the photocathode that is posi- 
tioned in front of the diode array (Heap et al. 1995). 

3. THE GHRS SPECTRA OF a CENTAURI A AND a 
CENTAURI B 

3.1. Data Analysis 

The nearby a Cen star system has been extensively 
studied in the past. In the ultraviolet, both a Cen A and B 
have been popular targets for the IUE. One of the more 
extensive analyses of IUE observations of a Cen A and B is 
that of Jordan et al. (1987), who used the IUE data to create 
simple one-dimensional models of the atmospheric struc- 
ture of the two stars. Ayres et al. (1995) and Hallam, Altner, 
& Endal (1991) have investigated the time variation of the 
UV line fluxes of a Cen A and B. Of particular interest is the 
substantial long-term variability noted by Ayres et al. (1995) 
for a Cen B, which may indicate a magnetic cycle similar to 
the 11 yr solar cycle and the activity cycles observed on 


TABLE 2 


GHRS Observations 


Star 

Grating 

Aperture 

Spectral Range 

(A) 

Spectral Resolution 
(km s -1 ) 

Date 

Start Time 
(UT) 

Exposure Time 

(s) 

Number of 
Readouts 

a Cen A 

G140M 

LSA 

1390-1417 

11 

1995 May 1 

14:38 

485 

2 


Ech-A 

SSA 

1212-1219 

3.5 

1995 May 1 

15:02 

2786 

12 


Ech-B 

SSA 

2792-2807 

3.5 

1995 May 1 

17:40 

485 

1 

a Cen B 

Ech-B 

SSA 

2792-2807 

3.5 

1995 May 5 

10:25 

700 

1 


G140M 

LSA 

1390-1417 

11 

1995 May 5 

11:44 

700 

2 


Ech-A 

SSA 

1212-1219 

3.5 

1995 May 5 

12:04 

3231 

12 

p Cet 

G160M 

LSA 

1385-1421 

15 

1994 Jun 2 

22:58 

2585 

8 

G160M 

LSA 

1537-1573 

15 

1994 Jun 3 

00:36 

1293 

4 

31 Com 

G160M 

LSA 

1385-1421 

15 

1994 Nov 7 

22:39 

3878 

12 


G160M 

LSA 

1537-1573 

15 

1994 Nov 8 

01:30 

1293 

4 

P Dra 

G160M 

LSA 

1532-1568 

20 

1992 Apr 23 

20:55 

889 

3 


G160M 

LSA 

1532-1568 

20 

1993 Feb 20 

21:15 

889 

3 


G160M 

LSA 

1383-1419 

15 

1995 Apr 29 

22:23 

1777 

6 


G160M 

LSA 

1532-1568 

15 

1995 Apr 30 

00:09 

889 

3 

P Gem 

G160M 

LSA 

1532-1568 

15 

1995 Sep 20 

02:08 

1293 

4 

G160M 

LSA 

1384-1420 

15 

1995 Sep 20 

03:29 

2047 

4 
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Fig. 1.- — GHRS spectra of the Si rv spectral regions of a Cen A and a Cen B, with line identifications 


many other solar-like stars using optical observations of the 
Ca ii H and K lines (Baliunas et al. 1995). 

Our GHRS a Cen A and B data sets unfortunately are 
very small, although the quality of the data is much 
superior to the IUE data. The G140M spectra of the 1390- 
1417 A spectral regions of a Cen A and a Cen B are shown 
in Figure 1. Besides the Si iv and O iv] lines, which were the 
main objectives for these observations, we also have identi- 
fied several weak chromospheric lines using the solar spec- 
tral atlas of Sandlin et al. (1986). We measured the 
velocities, widths (FWHM), and fluxes of these lines using 
single-Gaussian fits. The best fits were determined by % 2 
minimization (see, e.g., Bevington & Robinson 1992). The 
resulting parameters are listed in Tables 3 and 4. The 1 o 
random errors in the tables were estimated using Monte 
Carlo techniques. For all of the fits in this paper, we correct 
for instrumental broadening by convolving each proposed 


fit with the instrumental profile, which we assume is a 
Gaussian with a FWHM of 4.4 pixels (Gilliland 1994). 

Besides the lines identified in the 1390-1417 A spectral 
region, we also list the parameters of the H i, O v], and Mg 
n lines observed in echelle spectra analyzed by Linsky & 
Wood (1996). The self-reversed H i Lyman-a and Mg n h 
and k lines have high opacities and thus are not close to 
Gaussian in shape. We estimated velocities for these lines 
according to the centroid of the line wings, and we esti- 
mated fluxes using corrections for interstellar absorption 
from Linsky & Wood (1996). The O v] A1218 lines are 
shown explicitly in Figure 2. 

Laboratory wavelengths in Tables 3 and 4 are from Kelly 
(1987) and Sandlin et al. (1986), except for the O v] 11218 
wavelength, which is from Brown (1980). Based on the 
orbital parameters of Heintz (1982), we compute radial 
velocities of —23.4 and —20.3 km s -1 for a Cen A and a 


TABLE 3 

Line Profile Parameters for a Cen A 



Aest 

•^meas 

v a 

/ 

FWHM 


Ion 

(A) 

(A) 

(km s _1 ) 

(IQ" 13 ) 

(km s -1 ) 


Single-Gaussian Fits 

Hi b 

1215.670 

1215.571 

-1 + 1 

1200 + 200 



O v] 

1218.344 

1218.259 

2 + 4 

3.1 + 0.4 

65 + 9 

0.577 

Fe ii 

1392.817 

1392.707 

0 + 3 

0.12 + 0.03 

18 + 6 

0.957 

Si iv 

1393.755 

1393.673 

6 + 1 

9.5 + 0.1 

59 ± 1 

2.198 

O iv] 

1399.774 

1399.657 

-2 + 5 

0.21 + 0.04 

53 ±12 

1.109 

O iv] 

1401.156 

1401.082 

8 + 2 

0.98 + 0.06 

55 + 4 

0.858 

Si rv 

1402.770 

1402.681 

4+1 

4.7 + 0.1 

56 ± 1 

1.179 

O iv] 

1404.806 

1404.707 

2 ± 3 

0.42 + 0.05 

54 + 8 

1.036 

Ni n 

1411.071 

1410.955 

-1 + 3 

0.24 + 0.04 

38 + 7 

0.992 

Ni 

1411.949 

1411.848 

2 ± 3 

0.15 + 0.03 

27 + 7 

0.992 

Fe ii + Ni n + S i 


1412.745 


0.25 + 0.03 

26 + 5 

0.992 

Ni n 

1415.728 

1415.613 

-1 + 3 

0.14 + 0.03 

30 + 7 

1.040 

S rv]? 

1416.930 

1416.799 

-4 + 5 

0.16 + 0.03 

35 + 8 

1.040 

Mg n b 

2795.528 

2795.310 

0+1 

2250 + 80 



Mg n b 

2802.705 

2802.486 

0 + 1 

1730 + 80 



Two-Gaussian Fits 

Si rv 

1393.755 

1393.678 

7 + 1 

7.4 + 0.5 

47 ±2 

1.252 


1393.755 

1393.639 

— 2 ± 3 

2.5 + 0.4 

109 ± 10 

1.252 


a The line velocity relative to the star, assuming a stellar radial velocity of — 23.4 km s \ 
b The fluxes and velocities of the Lya and Mg n lines were measured directly rather than by a Gaussian fit. 
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TABLE 4 


Line Profile Parameters for a Cen B 



Aest 

X 

' w meas 

V* 

/ 

FWHM 


Ion 

(A) 

(A) 

(km s -1 ) 

(10- U) 

(km s" 1 ) 

xl 

Single-Gaussian Fits 

Hi b 

1215.670 

1215.588 

0 + 1 

1800 + 400 



O v] 

1218.344 

1218.275 

3 + 1 

1.8 + 0.1 

40 + 3 

0.813 

Si rv 

1393.755 

1393.710 

11 ± 1 

6.4 + 0.1 

50+ 1 

2.253 

O iv] 

1399.774 

1399.719 

9 + 4 

0.10 + 0.02 

40 + 9 

0.721 

O iv] 

1401.156 

1401.121 

13 + 1 

0.47 + 0.03 

42 + 3 

0.792 

Si iv 

1402.770 

1402.716 

9+1 

3.09 + 0.07 

48 + 1 

1.505 

O iv] ( + S iv]) 

1404.806 

1404.738 

6 + 2 

0.21 + 0.02 

38 + 6 

0.709 

S iv] 

1406.060 

1405.990 

5 + 6 

0.09 + 0.02 

45 + 11 

0.974 

O iv] 

1407.386 

1407.323 

7 + 6 

0.07 + 0.02 

43 + 11 

1.066 

Ni 

1411.949 

1411.884 

6 + 5 

0.04 + 0.02 

14 + 9 

0.884 

Fe n + Ni n + S i 


1412.764 


0.05 + 0.02 

17 + 8 

0.884 

Mg n b 

2795.528 

2795.339 

0+1 

1640 + 60 



Mg n b 

2802.705 

2802.515 

0+ 1 

1210 + 60 



Two-Gaussian Fits 

Si iv 

1393.755 

1393.711 

11 + 1 

4.7 + 0.3 

39 + 2 

0.953 


1393.755 

1393.686 

5 + 2 

1.8 + 0.3 

86 + 6 

0.953 

Si iv 

1402.770 

1402.722 

11 + 1 

2.3 + 0.2 

39 + 2 

0.953 


1402.770 

1402.697 

5 + 2 

0.9 + 0.2 

86 + 6 

0.953 


a The line velocity relative to the star, assuming a stellar radial velocity of — 20.3 km s 1 . 
b The fluxes and velocities of the Lya and Mg n lines were measured directly rather than by a Gaussian fit. 


Cen B at the time of the observations. Tables 3 and 4 list the 
line velocities in the stellar rest frame. 

The wavelengths of the oc Cen B spectra should be very 
accurate, because these observations were obtained imme- 
diately following wavelength calibration images. The 
Lyman-a and Mg ii h and k line wings appear to be well 



1217.5 1218.0 1218.5 1219.0 

Wavelength 

Fig. 2. — Single-Gaussian fits to the O v] 21218 lines of a Cen A and B 
observed with the Ech-A grating of the GHRS. The data have been 
smoothed slightly for the sake of clarity. The parameters of these fits are 
listed in Tables 3 and 4. 


centered on the rest frame of the star, which confirms the 
accuracy of the wavelength calibrations. The wavelength 
scale of the G140M spectrum of a Cen A is more uncertain, 
however, because the wavelength calibration image 
attempted prior to the observation was not completed suc- 
cessfully. We were therefore forced to use the default wave- 
length calibration for this spectrum, although we did do an 
offset correction as described by Soderblom, Sherbert, & 
Hulbert (1993). 

Initially, we found that the Fe ii, Ni n, and N i lines of a 
Cen A appeared to be blueshifted relative to the star by an 
average of 3 km s -1 . The blueshifts are unlikely to be 
genuine, because on the Sun and solar-like stars, weak 
chromospheric lines such as these tend to be at the stellar 
rest velocity or slightly redshifted (see, e.g., Wood et al. 
1996; Achour et al. 1995; Hassler, Rottman, & Orrall 1991). 
Therefore, we attribute the apparent blueshifts to an inaccu- 
rate wavelength calibration, and we compensate for them 
by increasing the wavelengths of the G140M spectrum by 
0.014 A (i.e., 3 km s -1 ). The velocities in Table 3 are the 
corrected values. We note that the H i and Mg ii lines of a 
Cen A appear to be centered on the stellar radial velocity, 
which suggests that there is no problem with the wave- 
length calibration of the echelle data. 

The last column of Tables 3 and 4 lists the xl values for 
the fits. In this paper, a fit will be considered acceptable if it 
is a 2 a fit, which means that the probability that random 
errors cannot account for the discrepancies between the fit 
and the data is <95.4%. (If the probability is <68.3%, it is 
a 1 cr fit, and so on.) The exact value of xl that separates the 
acceptable fits (i.e., 2 a fits) from unacceptable fits depends 
on the number of degrees of freedom, v, which for our pur- 
poses is basically the number of data points in the spectral 
region being fitted (see Bevington & Robinson 1992). This 
value will generally be about xl = 1-25 ± 0.10. 

Inspection of the xl values in Tables 3 and 4 reveals that 
the single-Gaussian fits to the stronger Si iv line of a Cen A 
and both Si iv lines of a Cen B are not acceptable. There- 
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fore, following the example of Linsky & Wood (1994), 
Linsky et al. (1995), and Wood et al. (1996), we fitted these 
lines with two Gaussians. In order to constrain the fits to 
the Si iv lines of a Cen B, the two lines were fitted simulta- 
neously, forcing the narrow and broad components of the 
two lines to have the same velocities and widths. For both 
the broad and narrow components, we also forced the Si iv 
2 1394/2 1403 flux ratio to be the optically thin value of 2. 
The two-Gaussian fits are shown in Figure 3, and the fit 
parameters are given in Tables 3 and 4. These results will be 
discussed further in § 5. 

3.2. Line Velocities 

Skylab and more recent rocket experiments have clearly 
demonstrated that the transition region lines of the Sun are 
redshifted relative to the solar photosphere and chromo- 
sphere (see, e.g., Doschek, Feldman, & Bohlin 1976; Hassler 
et al. 1991; Brekke 1993). The IUE has discovered similar 
line redshifts on many other stars, including a Cen A and B 
(Ayres, Jensen, & Engvold 1988). On the Sun, the observed 
redshifts are correlated with line formation temperature, 



Velocity (km s *) 


Fig. 3. — Two-Gaussian fits to the Si rv lines of a Cen A and B. The data 
are shown in histogram form. For each fit, the individual Gaussian com- 
ponents are displayed as thin solid lines, and the thick solid line is the 
convolution of the sum of the components with the instrumental profile. 
The parameters of these fits are listed in Tables 3 and 4. 


with a steady increase of line redshift up to about 10 km s“ 1 
at log T = 5.2. The behavior of the redshifts above 
log T = 5.2 is somewhat controversial. Based on the N v 
221239,1243 and O v] 21218 lines, which are formed at 
about log T — 5.3, there is evidence that the redshifts 
decrease dramatically above log T = 5.2. However, Brekke 
(1993) pointed out that this apparent decrease in redshift 
could be spurious due to uncertainties in the rest wave- 
lengths of the aforementioned N v and O v] lines. Indeed, 
the solar O v 21371 velocity is not consistent with the N v 
and O v] velocities, and suggests at most only a small 
decrease in redshift for temperatures above log T = 5.2. 
This issue is an important one, because the origin of the 
transition region line redshifts remains controversial, and 
the dependence of the redshifts on line formation tem- 
perature provides a potentially useful test of theoretical 
models that seek to explain the redshifts. 

The a Cen A and B line lists in Tables 3 and 4 may be 
short, but the data have unprecedented quality, and the 
lines span a broad range of temperature — from log T = 4.1 
for the Mg n lines to log T = 5.3 for the O v] 21218 line. 
Thus, we can at least crudely determine the variation of line 
redshifts with temperature. Figure 4 depicts the line red- 
shifts of a Cen A and B as a function of formation tem- 
perature. We plot the velocity of only the strongest O iv] 
line, and we restrict consideration to the measurements with 
errors less than +4 km s -1 . The line formation tem- 
peratures were estimated from the centroids of the contribu- 
tion functions (see, e.g., Brown et al. 1984; Jordan et al. 
1987) and the ionization equilibrium calculations of Arnaud 
& Rothenflug (1985). The line redshifts of the Sun and 
Procyon are also shown in Figure 4, based on data pre- 
sented by Achour et al. (1995) and Wood et al. (1996). 

Below temperatures of about log T = 5.2, the line red- 
shifts of the Sun, a Cen A, a Cen B, and Procyon are all very 
similar, although the Si rv and O rv] redshifts of a Cen B 
appear to be roughly 4 km s~ 1 larger than for the other 
three stars. Above log T = 5.2, the O v] 21218 lines of both 
a Cen A and B suggest drop-offs in redshift similar to that 
suggested by the solar data, when the same O v] rest wave- 
length is assumed. This is not the case for Procyon, however. 
The rest wavelengths of the N v and O v] lines, which 
suggest drop-offs in the line redshifts above log T — 5.2 for 
the Sun, a Cen A, and a Cen B, are 1238.821 A and 1218.344 
A respectively. Using these same rest wavelengths, we find 
no drop-off in line redshifts for Procyon, as previously noted 
by Wood et al. (1996). This provides strong evidence that 
the decrease in line redshift seen for the Sun, a Cen A, and a 
Cen B is real and not due to inaccuracies in the rest wave- 
lengths. The fact that the redshift decrease is seen for the 
Sun, a Cen A, and a Cen B, but is not seen for Procyon, a 
hotter and somewhat more active star, might represent an 
important insight into the origin of the effect. 

Achour et al. (1995) list a velocity of 9.0 km s -1 for the 
solar O v 21371 line. We did not plot this velocity in Figure 
4 because this line velocity is inconsistent with the N v 
21239 and O v] 21218 velocities. If the rest wavelengths of 
N v and O v] are accurate, as we have argued, then the rest 
wavelength of the O v 21371 line (1371.292 A) is suspect. 

3.3. T ransition Region Densities 

The O iv] lines seen in Figure 1 are very useful as density 
diagnostics (Cook et al. 1995). All four of the O rv] lines are 




No. 2, 1997 


EVALUATING POSSIBLE HEATING MECHANISMS 


751 



4 0 4.5 5.0 5.5 
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Fig. 4. — Line velocities of a Cen A, a Cen B, the Sun, and Procyon plotted as a function of line formation temperature. The data for the Sim and Procyon 
are from Achour et al. (1995) and Wood et al. (1996), respectively. The atomic species associated with each data point is indicated in the figure. In a number of 
instances, the velocities of two lines of the same atomic species are averaged, and only one data point is used to represent the lines in the figure. In these 
instances, the number of lines represented by the data point is indicated by a number in parentheses. 


detected in the a Cen B spectrum, while only three are 
clearly visible for a Cen A. For the latter, there appears to 
be a flux excess near 1407 A that probably is the O iv] 
1407.4 A line, but we cannot make an accurate measure- 
ment of this feature owing to the uneven continuum and 
possible line blends. However, we do have a flux measure- 
ment for the O iv] 1399.8 A line, which should have the 
same flux as the 21407.4 line for all densities (Cook et al. 
1995). In the midst of the O iv] multiplet falls the S iv] 
1406.1 A line, which is clearly visible in the a Cen B spec- 
trum. There is a feature at 1406 A in the a Cen A spectrum 
that probably is the S iv] line, but as was the case for O iv] 
21407.4, the uneven continuum and possible line blends 
prevent us from measuring the flux of the feature. Another 
complication is a S iv] line that is blended with O iv] 
21404.8. However, Cook et al. (1995) find that for relatively 
low densities, S iv] 21404.8 should have only 21% of the 
flux of S iv] 21406.1, which means that the O iv] 1404.8 A 


line accounts for at least 90% of the flux of the 1405 A 
feature in the spectra of both a Cen A and B. Nevertheless, 
for a Cen B we do reduce the O iv] 21404.8 flux reported in 
Table 4 by 10% (based on the measured S iv] 21406.1 flux) 
before proceeding with the analysis. 

There are three useful density-sensitive line ratios that 
can be computed from the O iv] lines: R 1 =/i 4 .oo//i 4 oi> 
^2 ~ /l40o//l405> an d ^3 ~/i 405 //i 40 i- For a Cen B, the 

average flux of the 21399.8 and 21407.4 lines is used for 
/ 1400 , since these two lines should have the same flux for all 
densities. From these three flux ratios and the O rv model of 
Cook et al. (1995), we compute three values for the electron 
density of each star. The flux ratios and corresponding den- 
sities are listed in Table 5. Note that for a Cen B, R 1 and R 2 
provide only upper limits for the density. The densities sug- 
gested by the three ratios are reasonably self-consistent, 
with R 3 providing the most accurate density measurement 
for both stars. Collectively, the results in Table 5 suggest 
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TABLE 5 

Electron Densities prom O iv] Line Ratios 


Star Line Ratio 2 log n e 


a Cen A R 1 = 0.21 ± 0.04 10.141?;^ 

R 2 = 0.55 ± 0.13 9.79±g;|i 

R 3 = 0.39 ± 0.06 9.6212:11 

a Cen B R 1 = 0.18 + 0.05 9.59±^ 80 

r 2 = 0.45 ± 0.15 9.56l£ 31 

R 3 = 0.40 ± 0.07 9.5612:4? 


2 R t ~/i40oZ^1401> ^2 — /l400/7l405> an< ^ ^3 ~ fl40s/fl*0V 

electron densities at log T = 5.14 of log n e = 9.65 ± 0.20 
and log n e = 9.50 ± 0.30 for a Cen A and a Cen B, respec- 
tively, and pressures (in units of cm -3 K) of log P e = 14.79 
± 0.20 and log P e = 14.64 ± 0.30. 

For comparison, the atmospheric models of a Cen A and 
B created by Jordan et al. (1987) ffrom IUE and Einstein 
data have transition region pressures of log P e = 14.80 and 
log P e = 15.28, respectively. Their a Cen A value agrees 
very well with our results, but their a Cen B pressure is 
about a factor of 4 higher than what we derive. Based on 
observations of a Cen made by the Extreme Ultraviolet 
Explorer (. EUVE ), Mewe et al. (1995) used line ratios of 
Fe x-Fe xiv lines to derive a best estimate for the electron 
density at about log T = 6.15 of log n e « 8.7 ± 0.5, which 
corresponds to a coronal pressure of log P e ~ 14.85 + 0.50. 
Since EUVE cannot resolve the two stars of the a Cen 
system, this pressure represents an average for the two stars. 
The inferred coronal pressure is consistent with the tran- 
sition region pressures we obtained from the O rv] lines. In 
contrast, based on similar comparisons of pressures derived 
from EUVE and O rv] data, Procyon and Capella appear 
to have coronal pressures that are significantly higher than 
their transition region pressures (Wood et al. 1996; Schmitt, 
Haisch, & Drake 1994; Linsky et al. 1995; Dupree et al. 
1993). This discrepancy may indicate that the line ratio 
diagnostics for the coronal plasma are unreliable, or it may 
indicate that on Procyon and Capella, the transition region 


line emission is to a large extent arising from regions physi- 
cally unconnected with those responsible for the coronal 
X-ray and EUV emission. Apparently, this does not have to 
be the case for a Cen A and B. 

4. THE C IV AND Si IV LINES OF 31 COMAE, p CETI, p 
DRACONIS, AND p GEMINORUM 

4.1. Data Analysis 

For 31 Com, P Cet, P Dra, and p Gem, many spectral 
regions were observed with the GHRS, but in this paper we 
will confine our attention to the C iv and Si iv lines. We 
fitted these lines with Gaussians, as we did for the Si iv lines 
of a Cen A and B in § 3.1. In addition to correcting for 
instrumental broadening, we also compensated for rota- 
tional broadening for 31 Com and P Dra, because both 
stars have v sin i values that are large compared with the 
spectral resolutions of the observations (see Tables 1 and 2). 
The parameters of the single-Gaussian fits are given in 
Table 6. The velocities in the table are in the stellar rest 
frames, based on the radial velocities listed in Table 1. For P 
Gem, the single-Gaussian fits, illustrated in Figure 5, are 
acceptable within the 2 a criterion. However, that is not the 
case for the single-Gaussian fits to the Si iv and C rv lines of 
the other three stars. Thus, the Si iv and C rv lines of 31 
Com, p Cet, and ft Dra were modeled with two Gaussians, 
as depicted in Figures 6, 7, and 8. The parameters of the fits 
are given in Table 7. 

In order to reduce the number of free parameters, we 
fitted the two C rv lines simultaneously, which makes the 
parameters of the weaker line dependent on those of the 
stronger line. For 31 Com, this is particularly useful because 
the two C iv lines are partially blended. By fitting the lines 
simultaneously and assuming the two lines have essentially 
the same profiles, the uncontaminated red wing of the 
weaker line constrains the shape of the red wing of 
the stronger line, and the uncontaminated blue wing of the 
stronger line constrains the shape of the blue wing of 
the weaker line. 

For p Ceti, fitting the lines simultaneously is also very 
useful because the two C iv lines are blended with at least 


TABLE 6 


Single-Gaussian Fit Parameters for 31 Com, ft Cet, p Dra, and fi Gem 


Star 

Ion 

Aest 

(A) 

Aneas 

(A) 

V* 

(km s" 1 ) 

/ 

(10~ 13 ) 

FWHM 
(km s ~ *) 

Xv 

jS Cet 

Si rv 

1393.755 

1393.874 

13 + 1 

6.55 + 0.06 

103 + 1 

5.206 

p Cet 

Si rv 

1402.770 

1402.880 

10+1 

3.29 + 0.04 

100+1 

2.414 

ft Cet 

Civ 

1548.202 

1548.347 

15 + 1 

5.50 + 0.09 

106 + 2 

1.892 

p Cet 

Civ 

1550.774 

1550.912 

14 + 1 

2.82 ± 0.07 

98 + 2 

1.416 

31 Com 

Si iv 

1393.755 

1393.779 

6 + 1 

4.70 + 0.04 

196 + 2 

1.829 

31 Com 

Si iv 

1402.770 

1402.764 

0+1 

2.59 + 0.03 

191 + 3 

1.309 

31 Com 

C iv 

1548.202 

1548.213 

3 + 1 

12.5 + 0.1 

240 + 3 

1.800 

31 Com 

Civ 

1550.774 

1550.815 

9 + 2 

6.4 + 0.1 

207 + 4 

1.800 

P Dra (1995) 

Si rv 

1393.755 

1393.779 

25 + 1 

7.68 + 0.08 

158 + 2 

1.568 

p Dra (1995) 

Si rv 

1402.770 

1402.734 

13 + 1 

4.51 + 0.07 

153 ±2 

1.204 

P Dra (1992) 

Civ 

1548.202 

1548.125 

5 + 1 

16.4 + 0.2 

195 + 2 

1.898 

p Dra (1992) 

Crv 

1550.774 

1550.753 

16 + 1 

9.0 + 0.1 

152 + 3 

1.898 

p Dra (1993) 

C iv 

1548.202 

1548.146 

9 ± 1 

16.7 + 0.2 

184 + 2 

1.355 

P Dra (1993) 

Crv 

1550.774 

1550.775 

20+ 1 

8.8 + 0.1 

141 + 3 

1.355 

p Dra (1995) 

C iv 

1548.202 

1548.187 

17+ 1 

13.4 + 0.2 

183 + 2 

1.521 

P Dra (1995) 

Civ 

1550.774 

1550.788 

23 + 1 

7.8 + 0.1 

148 + 3 

1.521 

p Gem 

Si rv 

1393.755 

1393.816 

10 + 1 

1.44 + 0.04 

82 + 2 

1.250 

P Gem 

Si iv 

1402.770 

1402.841 

12 + 1 

0.81 + 0.03 

80 + 3 

1.054 

p Gem 

C rv 

1548.202 

1548.300 

16 + 1 

1.63 + 0.07 

90 + 3 

1.133 

P Gem 

C iv 

1550.774 

1550.859 

13 ± 1 

0.82 + 0.03 

74 + 3 

1.219 


The line velocity relative to the star. 


No. 2, 1997 


EVALUATING POSSIBLE HEATING MECHANISMS 


753 




Wavelength 


F IG . 5. — Single-Gaussian fits to the Si rv and C rv lines of /? Gem. The 
parameters of these fits are given in Table 6. The Gaussian fit is shown 
before ( narrow solid line ) and after ( thick solid line ) convolution with the 
instrumental profile. The emission lines shown are (a) the Si rv A1394 line, 
(b) the Si rv A1403 line, and (c) the C rv XX 1548,1551 lines. A Si i line in the 
red wing of the C rv X 1551 line is also fitted with a Gaussian. 

three Si i emission features, which are fitted simultaneously 
with the C iv lines (see Fig. 6c). The Si i features, which were 
identified using the solar spectrum (Sandlin et al. 1986), 
obscure the far blue wing of the stronger line and the red 
wing of the weaker line. Fitting the C iv lines simulta- 
neously and assuming the two lines have essentially the 
same profile produces a more accurate separation of 
the C rv lines from the Si i lines. (Note that the Si i line at 
1551.3 A is also seen and corrected for in the /? Gem spec- 
trum shown in Fig. 5c.) 

When the C iv lines of HR 1099 and the Si iv lines of a 
Cen B were fitted simultaneously by Wood et al. (1996) and 
in § 3.1, respectively, the number of free parameters was 
reduced by forcing the widths and velocities of the narrow 
and broad components to be the same for both lines and by 
forcing the 21548/21551 and 21394/21403 flux ratios to be 
the optically thin value of 2. In other words, the two line 
profiles were assumed to be identical with the stronger line 
having exactly twice the flux of the weaker line. Unfor- 
tunately, we cannot assume this for the C iv lines of 31 Com 
and (1 Dra because the single-Gaussian fit parameters 



Wavelength 


Fig. 6. — Two-Gaussian fits to the Si rv and C rv fines of /? Cet. The 
parameters of these fits are given in Table 7. For each fit, the Gaussian 
components are plotted as thin solid lines, and the convolution of the sum 
of the components with the instrumental profile is shown as a thick solid 
line. The e mi ssion lines shown are (a) the Si rv 11394 line, (b) the Si rv 11403 
line, and (c) the C rv 11 1548,1551 fines. Two Si i features in the blue wing of 
the C rv 11548 fine and an additional Si i fine in the red wing of C rv 11551 
are also fit with Gaussians. 


clearly demonstrate that the line profiles are not identical. 
The stronger C iv line is broader and has a smaller redshift 
than the weaker line, and for ft Dra the 21548/21551 flux 
ratio is significantly less than 2. Linsky et al. (1995) found 
that the C iv lines of Capella have similar properties. The 
only way we can explain these properties of the C iv lines is 
to appeal to opacity effects and/or unidentified line blends. 
The implications of this conclusion will be discussed further 
in §§ 4.2 and 5. 

Thus, for the two-Gaussian fits to the C rv lines of 31 
Com and ft Dra, we allowed the velocities of the 21551 
narrow and broad components to be different from the 
velocities of the 21548 narrow and broad components, but 
we constrained the fit by forcing the velocity difference 
between the two broad components to be the same as the 
velocity difference between the narrow components. Simi- 
larly, we allowed the widths of the 21551 narrow and broad 
components to be different from the widths of the 21548 
narrow and broad components, but we forced the 
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Fig. 7. — Two-Gaussian fits to the Si iv and C rv lines of 31 Com. The 
parameters of these fits are given in Table 7. For each fit, the Gaussian 
components are plotted as thin solid fines, and the convolution of the sum 
of the components with the instrumental and rotational broadening pro- 
files is shown as a thick solid line. The emission lines shown are (a) the Si iv 
At 394 line, (b) the Si iv A 1403 fine, and (c) the C rv AA 1548, 1551 fines. The 
O rv] A1401 line is visible in the blue wing of the Si iv A1403 fine, and it has 
also been fitted with a Gaussian. 


FWHM A155l /FWHM A1548 ratio to be the same for both 
the narrow and broad components. Also, we required the 
21548/21551 flux ratio to be the same for both the broad 
and narrow components, although we allowed this ratio to 
be different from 2. Although there is no physical basis for 
these assumptions, they allowed us to lower the number of 
free parameters from 12 to 9, and the resulting fits are still 
excellent. 

Partially because of the Si i blends, it was not imme- 
diately apparent whether the C rv lines of /? Cet have the 
same profile or if the two lines are significantly different, like 
the C iv lines of Capella, 31 Com, and /? Dra. Thus, we fitted 
the C iv lines of /? Cet in the same manner as the C iv lines 
of 31 Com and fi Dra. The resulting parameters in Table 7 
demonstrate that for /? Cet, the two C iv lines are essentially 
identical and have a flux ratio very close to 2. 

The two Si iv lines of j B Cet and ft Dra were fitted simulta- 
neously in the same manner as the two C rv lines. For 31 
Com, however, the 21403 line is blended with the O iv] lines 




1546 1548 1550 1552 

Wavelength 


Fig. 8. — Two-Gaussian fits to the Si rv and C iv lines of /? Dra. The 
parameters of these fits are given in Table 7. For each fit, the Gaussian 
components are plotted as thin solid lines, and the convolution of the sum 
of the components with the instrumental and rotational broadening pro- 
files is shown as a thick solid fine. The emission fines shown are (a) the Si iv 
A1394 line, ( b ) the Si rv A1403 line, and (c) the C iv AA 1548, 1551 fines 
observed in three different years, where the dashed lines indicate the zero 
flux levels for each spectrum. The O iv] A1401 fine is visible in the blue 
wing of the Si rv A1403 fine, and it has also been fitted with a Gaussian. 

(especially the 21401 line; see Fig. lb) to the extent that we 
decided to fit the uncontaminated 21394 line independently 
of the 21403 line. The 21403 line was then fitted simulta- 
neously with the blended O iv] lines, and the narrow and 
broad component parameters of the 21403 line were con- 
strained by the previously measured 21394 parameters in a 
manner analogous to that described above for the C rv 
lines. 

4.2. The Variable C iv Lines of Draconis 

HST first observed /? Dra on 1992 April 23. Two of the 
planned GHRS observations were G160M exposures of the 
1383-1419 A and 1532-1568 A spectral intervals, the former 
containing the Si rv lines and the latter containing the C iv 
lines. However, the Si rv observation failed owing to a mal- 
function of the GHRS grating carrousel, although the other 
observations were completed successfully. In order to com- 
plete the data set, /I Dra was observed again about a year 
later on 1993 February 20. In addition to the missing Si iv 
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TABLE 7 


Two-Gauss ian Fit Parameters for 31 Com, fi Ckt, and p Dra 




^rest 

X 

'‘'me as 

V * 

/ 

FWHM 


Star 

Ion 

(A) 

(A) 

(km s" 1 ) 

(10" 13 ) 

(km s~ l ) 

Xv 

P Cet 

Si iv 

1393.755 

1393.859 

9+1 

4.3 + 0.2 

74 + 2 

1.271 


1393.755 

1393.920 

22 + 2 

2.6 + 0.2 

185 + 6 

1.271 

p Cet 

Si iv 

1402.770 

1402.871 

9 + 1 

2.18 + 0.08 

73 + 2 

1.271 


1402.770 

1402.932 

22 + 2 

1.32 + 0.08 

181 + 6 

1.271 

p Cet 

Civ 

1548.202 

1548.332 

12+1 

3.4 + 0.2 

74 + 2 

1.047 


1548.202 

1548.395 

24 + 3 

2.4 + 0.2 

183 + 9 

1.047 

P Cet 

Civ 

1550.774 

1550.901 

11 + 1 

1.8 + 0.1 

70 + 2 

1.047 


1550.774 

1550.963 

23 + 3 

1.2 + 0.1 

175 + 11 

1.047 

31 Com 

Si rv 

1393.755 

1393.789 

8 + 2 

3.0 + 0.3 

144 + 7 

1.040 



1393.755 

1393.742 

-2 + 5 

1.9 + 0.3 

321 + 21 

1.040 

31 Com 

Si iv 

1402.770 

1402.773 

2 + 2 

1.6 + 0.1 

141 + 8 

1.006 



1402.770 

1402.726 

-8 + 5 

1.1 + 0.1 

314 + 22 

1.006 

31 Com 

Civ 

1548.202 

1548.232 

7 + 2 

5.0 + 06 

131 ±9 

1.153 



1548.202 

1548.188 

-2 + 3 

8.0 + 0.5 

325 ± 13 

1.153 

31 Com 

Civ 

1550.774 

1550.849 

16 + 2 

2.5 + 03 

109 + 8 

1.153 



1550.774 

1550.805 

7 + 4 

4.0 + 0.3 

270 + 11 

1.153 

p Dra (1995) 

Si iv 

1393.755 

1393.782 

26 + 1 

6.1 ± 0.4 

137 + 4 

1.045 


1393.755 

1393.765 

22 + 6 

1.8 + 0.4 

282 + 27 

1.045 

P Dra (1995) 

Si iv 

1402.770 

1402.736 

13 + 1 

3.6 + 02 

134 + 4 

1.045 


1402.770 

1402.719 

9 + 6 

1.0 + 0.2 

277 + 26 

1.045 

P Dra (1992) 

Civ 

1548.202 

1548.128 

6 + 1 

11.4 + 0.8 

150 + 6 

1.240 


1548.202 

1548.132 

7 + 6 

5.7 + 0.8 

353 + 30 

1.240 

P Dra (1992) 

Civ 

1550.774 

1550.752 

16 + 2 

6.2 + 0.5 

115 + 5 

1.240 


1550.774 

1550.756 

17 + 6 

3.1 + 0.4 

272 + 23 

1.240 

P Dra (1993) 

Civ 

1548.202 

1548.154 

11 + 1 

12.1 + 1.3 

153 ±7 

1.149 


1548.202 

1548.113 

3 + 6 

4.9 + 1.2 

291 + 27 

1.149 

P Dra (1993) 

Civ 

1550.774 

1550.780 

21 + 2 

6.4 + 0.7 

117 + 5 

1.149 


1550.774 

1550.738 

13 + 6 

2.6 + 0.6 

223 + 20 

1.149 

P Dra (1995) 

C rv 

1548.202 

1548.178 

16 + 2 

7.8 + 0.8 

134 + 7 

1.136 


1548.202 

1548.221 

24 + 4 

6.2 + 0.8 

277 + 18 

1.136 

P Dra (1995) 

C rv 

1550.774 

1550.770 

20 + 2 

4.5 + 0.5 

109 + 5 

1.136 


1550.774 

1550.814 

28 + 4 

3.6 + 05 

224 + 14 

1.136 


a The line velocity relative to the star. 


spectrum, another C rv spectrum also was planned. 
Although the C rv spectrum was obtained successfully, the 
Si iv observation failed again due to yet another carrousel 
malfunction. On 1995 April 29, /? Dra was acquired for a 
third time by HST for observations of Si rv and C rv, and 
finally both spectra were obtained without incident. One 
positive result of the repeated failure of the Si rv observation 
is that GHRS spectra of the C rv lines were acquired in 
three different years, which allowed us to check for fine 
profile variability. 

Figure 9 illustrates the smoothed C rv profiles for all 
three years, and Figure 8c shows the two-Gaussian fits to 



Wavelength 


F IG . 9. — Fourier-smoothed spectra of the C rv lines of ff Dra observed 
in three different years. Line profile variability is apparent. 


these lines. Line profile variability is evident in both Figure 
9 and in the C rv parameters listed in Tables 6 and 7. The 
C rv fluxes appear to be constant between 1992 and 1993, 
but the 21548 and 21551 fluxes dropped by 18% and 10%, 
respectively, between 1993 and 1995. The effects of inaccu- 
racies in the GHRS flux calibration must be considered 
here, because the 1992 and 1993 data are pre-COSTAR 
observations, while the 1995 spectrum was taken after the 
installation of COSTAR, which affected the instrumental 
sensitivities. However, flux calibration inaccuracies cannot 
explain why the flux decrease is different for the two C iv 
lines, so we believe the flux decrease is due to true stellar 
variability. Since the variability is more pronounced in the 
more optically thick C rv 21548 line, opacity effects prob- 
ably play a significant role in altering the C rv profiles. 

In addition to the flux change, the 1995 C rv lines also are 
more redshifted than they were in 1992 and 1993, and once 
again the change is larger for the 21548 line. An increase in 
line redshift also is evident between 1992 and 1993. Quanti- 
tatively, the measured 21548 (21551) line velocities based on 
the single-Gaussian fits are 5 + 1, 9 + 1, and 17 + 1 
(16 ± 1, 20 + 1, and 23 ± 1) for the 1992, 1993, and 1995 
data, respectively. The velocity variability is too large to be 
due to inaccuracies in the wavelength calibration. 

The obvious velocity and flux variability, combined with 
more subtle line profile changes, naturally have affected the 
parameters of the two-Gaussian fits. For example, the wings 
of the 1992 C rv profiles appear to be slightly more extended 
than is the case for the other two years, resulting in wider 
broad components. Also, the upper halves of the 1995 C rv 
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profiles are narrower than is the case for the other two 
years, whch results in thinner narrow components. The 
broad component contribution to the C iv lines also is 
significantly larger for the 1995 profiles than it is for the 
1992 and 1993 profiles. 

As discussed in § 1, we currently are working under the 
hypothesis that the broad components are associated with 
microflaring. Thus, the larger broad component contribu- 
tion to the 1995 C iv profiles appears to suggest a greater 
prevalence of microflaring in 1995 than in 1992 and 1993. 
However, we mentioned above that opacity effects probably 
play some role in the line variability. Opacity effects and/or 
unidentified line blends also appear to be necessary in order 
to explain the general velocity, width, and flux discrepancies 
between the two C iv lines of /? Dra, Capella, and 31 Com 
(see above). To some extent, these effects compromise the 
ability of the two-Gaussian fits to accurately quantify the 
narrow and broad component contributions to the line 
emission. We will return to this issue in § 5. 

With only three /? Dra C rv observations spaced so far 
apart in time, we cannot determine whether the detected 
variability is associated with different regions of the star 
rotating into and out of view, or whether the variability is 
more global, perhaps like that associated with the 11 yr 
activity cycle of the Sun. Furthermore, because the profile 
changes primarily are on the blue sides of the C iv profiles, 
we cannot rule out a warm stellar wind as the source of the 
variability. The radius and v sin i values in Table 1 suggest 
that the rotational period of ft Dra is roughly 6 months, 
which is less than the time between the observations. Thus, 
rotational modulation is a plausible explanation. 

5. COMPARING STELLAR TRANSITION REGION LINE 
PROFILES 

5.1. Broad Components and Microflaring 

The primary goal of this paper is to explore the variations 
of stellar transition region line profiles with activity and 
other stellar properties such as surface gravity. The tran- 
sition region lines for ten of the stars in our sample are 
displayed in Figure 10 along with the Gaussian fits to the 
profiles. For each star, the line displayed is either Si iv 
A1394 or C rv A1548, depending on which is the highest 
quality transition region line profile available. For most of 
the stars, two Gaussians were required to fit the line pro- 
files. For Capella, an additional third Gaussian was added 
to represent the emission from the G8 star of the binary 
system. As discussed in § 4.1, a second Gaussian was not 
required to model the Si iv line of p Gem accurately, so a 
single-Gaussian fit is displayed. 

For Procyon, a simple single-Gaussian fit is used despite 
the obvious presence of a blue wing excess. The blue wing 
excess is seen for all of Procyon’s transition region lines 
(Wood et al. 1996). However, except for Procyon’s N v 
A1239 line, there is no corresponding red wing excess that 
would indicate the presence of a broad component like that 
seen for most of the stellar Si iv and C iv lines in Figure 10. 
Attempts to fit two Gaussians to several of Procyon’s tran- 
sition region lines resulted in a reproduction of the single- 
Gaussian fit plus a very small Gaussian in the blue wing of 
the line, rather than a narrow component plus a broad 
component. Wood et al. (1996) therefore concluded that 
Procyon’s blue wing excesses have a different interpretation 
than the broad components observed in the transition 


region lines of the other stars. Since the single-Gaussian fits 
to Procyon’s transition region lines are acceptable 2 o fits 
when the blue wings are ignored, no two-Gaussian fits were 
presented, with the notable exception of N v A1239. 

We converted previously published X-ray fluxes for our 
sample stars to surface fluxes using the stellar radii and 
distances listed in Table 1. These X-ray surface fluxes are 
given in Table 8, together with the C iv surface fluxes (both 
lines combined). Both quantities will be considered to be 
measures of the stellar activity level in this paper. The stars 
in Table 8 are listed in order of decreasing X-ray surface 
flux, and the same order is used in Figure 10. Columns 
(4)-(9) of Table 8 list properties of the fits to the Si iv and 
C iv lines, and the last column indicates references for the 
quantities listed in the table (in addition to this paper). The 
ratio of the broad component flux to the total line flux is 
given in columns (4) and (5). The difference between the 
velocities of the broad and narrow components is given in 
columns (6) and (7). We converted the widths of the broad 
and narrow components to nonthermal velocities using the 
relation (in cgs units) 



= 3.08 x 10~ 21 




(1) 


where AX is the measured FWHM, T is the line formation 
temperature, m ; is the ion mass, and £ is the most probable 
nonthermal speed. The narrow and broad component non- 
thermal velocities listed in columns (8) and (9) of Table 8 
represent weighted averages of all four Si iv and C rv lines. 

Based on Figure 10, there appears to be a tendency for 
the transition region lines of active stars to have more non- 
Gaussian profiles with broader wings than the less active 
stars. A particularly useful exercise is to compare the Si iv 
profiles of P Cet and ft Gem (see also Figs. 5 a and 6a), 
because these two stars are very similar in all respects 
except activity level. The wings of the ft Cet Si rv fine clearly 
are more prominent than those of ft Gem. 

The profiles in Figure 10 are compared more directly in 
Figure 11, in which the C iv profile of AB Dor from Yilhu et 
al. (1996) is also included. The profiles shown in Figure 11 
are based on the two-Gaussian fits (except for ft Gem). The 
Procyon profile shown in Figure 1 1 is also based on a two- 
Gaussian fit, because for the purposes of this figure we 
wanted the aforementioned blue-wing excess to be evident. 
The stellar line profiles in Figure 11 have been normalized 
to have the same amplitude and FWHM (1.0 and 100 km 
s -i , respectively), and at the half-maximum location we 
force the lines to be centered on 0 km s _1 . A logarithmic 
scale for the flux is used to emphasize the line wings. The 
three stars with the largest X-ray surface flux (AU Mic, HR 
1099, and AB Dor) clearly have the broadest line wings, 
once again suggesting a correlation between activity and 
extended wing emission. 

Based on Figures 10 and 11, we have argued that a corre- 
lation exists between stellar activity and broad transition 
region line wings. The question we now address is, what is 
the best way to quantify the prominence of the line wings? 
We propose two methods. The first is to measure the width 
of a profile at 5% of the maximum amplitude and then to 
divide this width by the FWHM. In Figures 12a-12c, we 
plot this ratio against three traditional activity diagnostics 
for the chromosphere, transition region, and corona: the 
Mg ii, C iv, and X-ray surface fluxes, respectively. The Mg ii 
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Fig. 10. — Our best Gaussian fits to the Si iv or C rv lines of nearly all the stars in our sample. The emission line shown is either the Si iv A 1394 line or the 
C iv X 1548 fine, depending on which is the most free of blends and has the highest S/N. All of these fits have already been displayed before, either earlier in this 
paper or in previous papers (Linsky & Wood 1994; Linsky et al. 1995; Wood et al. 1996; Wood 1996), but here we show them all in one place to allow for 
easier comparison. With the exception of Procyon and /? Gem, two Gaussians are required to fit the data adequately. For Capella, a third Gaussian 
representing emission from the G8 III star is also included. The stars are placed in order of decreasing X-ray surface flux, reading from AU Mic to Procyon 
and then down from ft Cet to /? Gem. 


fluxes were measured from GHRS data, except for the AB 
Dor value which is from 1UE data (Rucinski 1985). Basal 
fluxes from Rutten et al. (1991) have been subtracted from 
the Mg ii fluxes. The second method of quantifying the 
prominence of the line wings is to divide the broad com- 
ponent flux by the total line flux. These values, which are 
listed in columns (4) and (5) of Table 8 for both Si iv and 
C iv, are plotted versus the Mg n, C rv, and X-ray surface 
fluxes in Figures \2d-\2f. 

Figures 12c and 12 d clearly indicate a correlation 
between broad line wings and stellar activity. This corre- 
lation is less obvious in Figures 12 b and 12e and is not very 
apparent at all in Figures 12a and 12 d. The reason that the 
correlation is most apparent when F x is used as the activity 
diagnostic is probably that the F x values vary by about 5 
orders of magnitude, while the F MgII an( ^ ^civ values vary 
by only about 2 and 3 orders of magnitude, respectively. 


Thus, scatter in the data points produced by systematic 
errors can hide the correlation in the Mg n and C iv plots 
more easily than in the X-ray plots. Variations in line profile 
behavior with spectral type are probably the largest source 
of scatter in Figure 12. In determining how line profiles vary 
with activity, ideally one should observe stars with the same 
spectral type but different activity levels. Unfortunately, 
such a data set is not available at this time, so our sample 
necessarily contains cool stars with a large range of spectral 
types. Other causes of the scatter in Figure 12 include 
uncertainties in F MgII , F CIV , and F x due to both stellar 
variability and uncertainties in stellar radii and distances, 
line profile variability such as that seen for Dra, and 
uncertainties introduced into our analysis by opacity effects 
and line blending, as discussed in § 4.2. 

The upper limit representing Procyon in Figures 12d-12f 
appears to be inconsistent with the general trend implied by 
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TABLE 8 


Transition Region Line Profile Properties 


Star 

(1) 

] og F x 
(2) 

Pc IV 

(3) 

/bc/7ioi 

U NC 

^BC 

£nc 

(km s *) 
(8) 

£bc 

(km s x ) 

(9) 

References 

(10) 

Si rv 
(4) 

C iv 
(5) 

Si rv 
(6) 

C IV 

(7) 

AU Mic 

7.55 

5.28 

0.56 + 0.07 

0.43 + 0.03 

2 + 5 

-9 + 3 

15.4 + 1.5 

96.5 + 9.2 

1, 2, 3 

HR 1099 

7.40 

5.84 

0.55 + 0.14 

0.66 + 0.03 

(0) 

8 + 2 

48.5 + 4.7 

149.7 + 2.4 

4, 5 

AB Dor 

6.79 

5.49 


0.49 + 0.03 


6 + 4 

39.1 + 1.2 

200.3 + 4.8 

6, 7 

31 Com 

6.22 

5.57 

0.39 + 0.06 

0.62 + 0.04 

10 + 5 

9 + 4 

78.8 + 9.7 

179.1 + 17.0 

8 

Capella 

5.38 

5.33 

0.35 + 0.04 

0.52 + 0.05 

19 ±4 

-4 + 4 

78.3 + 5.8 

186.9 + 13.4 

3, 8, 9 

Procyon 

5.19 

4.65 

<0.20 

<0.20 



45.3 + 1.3 


5, 8 

P Cet 

5.02 

3.76 

0.38 + 0.03 

0.41 ± 0.03 

— 13 + 2 

— 12 ± 3 

42.7 + 1.8 

109.1 + 4.1 

8 

a Cen B 

4.80 

3.81 

0.28 ± 0.05 


6 + 2 


22.6 + 1.2 

51.3 + 3.6 

8 

Sun 

4.51 

3.98 





26.0 + 5.0 


10, 11, 12 

p Dra 

4.38 

4.59 

0.23 + 0.05 

0.35 ± 0.04 

4 ± 6 

-1 + 8 

78.3 + 6.7 

160.0 + 20.4 

8 

a Cen A 

4.06 

3.73 

0.25 + 0.04 


9 + 3 


27.6 + 1.2 

65.2 + 6.0 

8 

p Gem 

3.16 

2.76 

<0.20 

<6.20 



48.2 + 3.4 


8 


References. — Wood et al. 1994; (2) Linsky & Wood 1994; (3) Wood 1996; (4) Dempsey et al. 1993a; (5) Wood et al. 1996; (6) Hempelmann et al. 
1995; (7) Vilhuet al. 1996; (8) Ayres et al. 1995; (9) Linsky et al. 1995; (lO)Maggio et al. 1987; (11) Ayres, Marstad, & Linsky 1981 ; (12)Dere & Mason 
1993. 


the other data points. However, Procyon might not be as 
discrepant as these figures suggest. The broad component 
contributions plotted in Figures 12d-12f are only for the 
Si rv and C iv lines. No broad component was detected 
for Procyon in these lines, but a broad component was 
detected in the N v X 1239 line. The N v broad component 
was found to account for 42% ± 8% of the total line flux, so 
if a N v data point were plotted for Procyon, it would be 
consistent with the Si iv and C rv data points of the other 
stars. The N v line is formed at somewhat higher tem- 
peratures than the Si rv and C rv lines. Perhaps for Procyon, 
which is the hottest star in our sample, the broad com- 
ponent phenomenon only becomes evident at higher tem- 
peratures in the transition region. 

If microflaring is indeed the correct explanation for the 
broad components, then Figure 12 suggests that the micro- 



-300 -200 -100 0 100 200 300 

Velocity (km s l ) 


Fig. 11. — Comparison of the transition region line profiles of all the 
stars in our sample, corrected for instrumental and rotational broadening, 
based on the fits to the data shown in Fig. 10 (except for Procyon — see 
text). The C iv profile of AB Dor is also shown, based on the two-Gaussian 
fit of Vilhu et al. (1996). In order to compare the profile shapes properly, all 
the lines have been normalized to have the same amplitude (1.0) and 
FWHM (100 km s -1 ), and at the half-maxi m um location we force the lines 
to be centered at a velocity of 0 km s -1 . The profiles have been plotted on a 
log flux scale in order to emphasize the line wings. The shaded region 
indicates what a Gaussian profile would look like when plotted this way. 
In the figure, the stars are listed in order of decreasing X-ray surface flux. 
The more active stars clearly tend to have the broader line wings. 


flaring must be more prevalent on active stars than on inac- 
tive stars. It certainly would not be surprising if this were 
the case. In addition to the microflaring hypothesis, Vilhu et 
al. (1996) proposed another interpretation for the broad 
component of AB Dor. They suggested that the broad wing 
emission could originate in large “ slingshot prominences ” 
that extend several stellar radii above the surface of the star 
and rigidly rotate with the star. Because AB Dor is rotating 
very rapidly (u sin i = 90 km s -1 ), emission in these sling- 
shot prominences would be rotationally broadened in a 
way that could produce the extensive wings observed in 
C iv. This hypothesis is plausible because Ha observations 
suggest that prominences this large do indeed exist on AB 
Dor (Cameron & Robinson 1989). It is uncertain, however, 
whether such structures contain gas hot enough and abun- 
dant enough to account for the broad components of the 
C rv lines. 

Figure 11 shows that the C rv lines of AB Dor and AU 
Mic have very similar shapes. However, it is unlikely that 
the slingshot prominence interpretation is also applicable 
for AU Mic, because its rotation rate is only v sin i < 5.9 
km s -1 (see Table 1). In order to explain the broad C iv 
wings of AU Mic, rigidly rotating prominences would have 
to extend out to at least 15 stellar radii, which seems 
improbable. Also, slingshot prominences cannot explain 
why the transition region line wings of /? Cet are more 
extended than those of the similar, but much less active star 
fi Gem, because both are slow rotators (see Table 1). Thus, 
we believe that microflaring is a better interpretation of the 
broad C rv line wings, as it seems to be more broadly applic- 
able. 

Another attraction of the microflaring interpretation is 
the possible existence of a well-studied solar analog of the 
phenomenon, namely transition region explosive events 
(see, e.g., Dere, Bartoe, & Brueckner 1989). The GHRS 
observations of a Cen A provide an excellent opportunity to 
test the connection between the solar explosive events and 
the stellar broad components because a Cen A is nearly 
identical to the Sun in almost all respects, including age 
(about 5 x 10 9 yr) and spectral type (G2 V). The width of 
the broad component of a Cen A is 109 + 10 km s" 1 , which 
agrees very well with the line widths typically produced by 
the solar explosive events. Although the centroids of the 
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Fig. 12. — Two quantities measuring the prominence of the line wing emission are plotted vs. the Mg n, C IV, and X-ray surface fluxes. The first quantity 
[A2(0.05)/AA(0.5)] is the width of the line profile at 5% of the maximum amplitude divided by the FWHM, which is plotted in panels (a)-(c). The second 
quantity (f BC /f at ) is the flux of the broad component of the two-Gaussian fit divided by the total line flux, which is plotted in panels (d)-(f) for both the C iv 
(boxes) and Si rv ( diamonds ) lines. 


emission lines produced by individual solar explosive events 
vary greatly, Dere et al. (1989) found that blueshifted pro- 
files are more common than redshifted profiles, which may 
explain why the broad component of a Cen A is blueshifted 
relative to the narrow component by 9 ± 3 km s “ 1 (see 
Table 8). Inspection of columns (6) and (7) of Table 8 reveals 
a general tendency for the broad components to be blue- 
shifted relative to the narrow components. The only clear 
examples of redshifted broad components are the Si iv and 
C iv lines of Ceti, and possibly the C rv lines of AU Mic. 

One problem with the solar explosive event interpreta- 
tion of the broad components is that the broad component 
of a Cen A accounts for 25% + 4% of the total line flux, 
whereas explosive events appear to account for only about 
5% of the Sun’s total transition region line flux (Linsky & 
Wood 1994). If the Sun were plotted in Figures 12c and 12 d 
with/ BC // tot = 0.05, it would be inconsistent with the general 
trend suggested by the other data points. There are a 
number of possible explanations for this problem : 

1. The amount of microflaring in the solar transition 
region is anomalously low. 

2. The prevalence of the explosive events in the solar 
transition region has been underestimated, possibly owing 
to insufficient spatial and/or time resolution. 

3. Stellar transition region explosive events are not 
responsible for the observed broad wing emission seen in 
the Si iv and C iv lines. 


4. The line fluxes produced by the stellar explosive events 
are overestimated by the two-Gaussian fits to the disk- 
averaged line profiles. 

The first of these explanations does not seem very likely, 
because to our knowledge the Sun is not a particularly 
unusual star in any other respect. It is difficult for us to 
assess the plausibility of the second explanation, but Parker 
(1988) has previously proposed that the explosive events 
observed by HRTS are only the most energetic of the solar 
microflares and that most of the microflaring activity on the 
Sun is as yet undetectable. Hopefully, the recently launched 
Solar Heliospheric Observatory (SOHO) will shed some light 
on this issue. If the third explanation is true, that would cast 
serious doubt on the whole microflaring interpretation of 
the broad wings, as we are aware of no other solar pheno- 
menon that could be the solar analog of the hypothetical 
microflaring. The last explanation is quite possible, because 
even if microflares are the correct explanation for the broad 
wings of the stellar transition region lines, this does not 
mean that simple two-Gaussian fits accurately quantify the 
properties of these microflares. 

If either the third or fourth explanations is true, then the 
nonmicroflaring component of the line profile must contrib- 
ute more emission to the line wings than our Gaussian 
narrow components suggest. We believe that this is almost 
certainly true to some extent for the following reasons. First 
of all, Dere & Mason (1993) find that solar-quiet region 
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Si iv and C iv line profiles generally have excess emission in 
the line wings when compared with a Gaussian profile fitted 
to the emission lines. (However, this might conceivably be 
due to unresolved microflaring.) Furthermore, even if we 
disregard explosive events, the widths of Si iv and C iv 
profiles still vary significantly from place to place on the 
solar surface. Even if we assume that these profiles are 
Gaussian, the sum of Gaussians with different widths will 
still result in a profile with wings extended beyond those of 
a Gaussian profile. Thus, even in the absence of transition 
region explosive events, we expect that disk-averaged solar/ 
stellar line profiles will have some broad wing emission. 
Wood et al. (1996) demonstrated that a decent fit to the C iv 
lines of HR 1099 could be achieved without a broad com- 
ponent, if one assumes that a very broad distribution of line 
widths are being emitted from the stellar surface (see their 
Fig. 17). We note here, however, that models such as this 
cannot explain the observed line asymmetries, which are 
indicated in Table 8 by the velocity differences between the 
broad and narrow components, seen for HR 1099 and most 
of the other stars in our sample. 

For the particular case of a Cen A, we have performed 
another two-component fit to the data in which we assume 
that the narrow component is not simply Gaussian but is 
instead emitted from a transition region producing Gauss- 
ian line profiles with a distribution of widths. We assume a 
normal distribution of line widths, with the distribution 
width parameter, a, set to be 30% of the average of the 
distribution. This assumed distribution roughly mimics the 
solar case (Dere & Mason 1993). The resulting fit is depicted 
in Figure 13. In addition to the two components and their 
sum ( solid lines), we also show a Gaussian profile ( dotted 
line) with the same velocity and flux as the narrow com- 
ponent of the fit and with a width identical to the average of 
the narrow component’s distribution of widths. Comparing 
the dotted line with the narrow component reveals the 



Velocity (km s l ) 

Fig. 13. — Two-component fit to the Si iv A 1394 line of « Cen A, where 
the two components are plotted as thin solid lines and the convolution of 
their sum with the instrumental profile is shown as a thick solid line. The 
broad component is simply a Gaussian, but we assume that the narrow 
component is the sum of many Gaussians with a normal distribution of 
widths. The width parameter of the distribution, a, is assumed to be 30% of 
the average of the distribution. This assumption is meant to consider the 
possibility that not all nonmicroflaring atmospheric regions produce line 
profiles with exactly the same width. The resulting narrow component is 
compared with a Gaussian profile ( dotted line) that has the same velocity 
and flux as the narrow component, and has a width equal to the average of 
the narrow component’s distribution of widths. 


extent to which the assumed line width distribution has 
broadened the wings of the narrow component. As a conse- 
quence, the broad component contribution to the line 
profile (fuc/fot) has dropped to 0.15 + 0.06 from a previous 
value of 0.25 ± 0.04 for the fit in Figure 3 a. This contribu- 
tion still is larger than the 5% value believed to be appropri- 
ate for the Sun, but this exercise does demonstrate how a 
simple two-Gaussian fit could significantly overestimate the 
broad component’s flux contribution to the Si iv and C iv 
lines. 

5.2. Narrow Components and Nonthermal Velocities 
In Figure 14a, we plot the transition region nonthermal 
velocities, as measured by the widths of single-Gaussian fits 
to the Si iv and C iv lines, versus the stellar surface gravity. 
The surface gravities were computed from the masses and 
radii listed in Table 1. In Figure 14b, we plot the non- 
thermal velocities of the narrow and broad components 
versus surface gravity. The dashed lines in both panels of 
Figure 14 indicate the sound speed at a temperature of 
10 5 K. The broad component’s nonthermal velocities in all 
cases are supersonic. This leads us to speculate that the 
large widths of the broad components are not due to turbu- 
lence, per se, because under normal circumstances super- 
sonic turbulence would be degraded rapidly by shock 
dissipation (see, e.g., Shu 1992). 



Fig. 14. — (a) Transition region nonthermal velocities, as measured 
using single-Gaussian fits to the Si rv and C rv lines, plotted vs. stellar 
surface gravity (normalized to the solar value). ( b ) Transition region non- 
thermal velocities, as measured by two-Gaussian fits to the Si rv and C rv 
lines, plotted vs. surface gravity for both the narrow ( diamonds ) and broad 
(boxes) components. The dashed lines indicate the sound speed at 10 s K. 
Also displayed is a linear least-squares fit to the narrow component veloci- 
ties of the six stars with the highest surface gravities ( solid line). 
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One of the attractive aspects of the microflaring interpre- 
tation of the broad components is that it suggests an alter- 
native explanation for the large line widths. Flares and 
microflares produce broad emission lines owing to high- 
velocity jets of plasma ejected from the magnetic reconnec- 
tion site rather than owing to isotropic turbulence. We do 
not expect line profiles produced under these conditions to 
be Gaussian, and in fact most solar explosive events gener- 
ate C iv lines that are far from Gaussian. The sum of thou- 
sands of such profiles, from different locations on the stellar 
disk, might result in a profile that is close to Gaussian, but 
the width of this Gaussian will not be characteristic of the 
true turbulent velocity. 

Of the six low-gravity, evolved stars in our sample, three 
have narrow components with widths that suggest non- 
thermal velocities comparable to the sound speed, and three 
((} Dra, Capella, and 31 Com) have narrow component non- 
thermal velocities about 60% larger than the sound speed 
(see Fig. 14b). The latter are the same three stars for which 
we believe that the C iv and Si iv lines are broadened by 
opacity effects and/or line blending (see § 4). Thus, it is likely 
that the narrow-component nonthermal velocities have 
been significantly overestimated for these three stars. We 
believe that our data in Figure 14b are consistent with the 
narrow component having a constant nonthermal velocity, 
equal to the sound speed, for stars with surface gravities 
between log g/g Q = —3.0 and log g/g Q = —0.4. For stars 
with larger gravities, the nonthermal velocity decreases with 
increasing surface gravity. We have performed an analytical 
least-squares fit to quantify this correlation for log g/g Q > 
- 0.4 (see Fig. 14b), and we find that 

( \ — 0.68 ± 0.07 

£) ' (2) 

Because of the efficiency of shock dissipation in damping 
supersonic motions, it is not surprising that the £ NC versus g 
relation flattens as soon as £ NC reaches the sound speed. 

Note that this apparently tight correlation between £ NC 
and g is not evident in Figure 14a. The solar nonthermal 
velocity plotted in Figure 14 is the result of spatially resolv- 
ed C iv and Si iv observations of solar quiet and active 
regions (Dere & Mason 1993). In Figure 14a the solar non- 
thermal velocity is lower than the stellar velocities mea- 
sured from the disk-averaged line profiles, but in Figure 14b 
the solar velocity agrees very well with the correlation 
between £ NC and g/g Q . The existence of the Figure 14b 
correlation, and the agreement of this correlation with the 
solar value, represents a possible validation of our use of 
two-Gaussian fits to quantify transition region line profiles. 

The significance of the narrow-component velocities is 
that they represent a substantial amount of kinetic energy, 
possibly enough to heat the entire transition region. Are 
these nonthermal velocities due to isotropic turbulence, or 
are they caused by the passage of waves through the tran- 
sition region? The flux in a propagating wave can be 
written as (see, e.g., Brown & Jordan 1981) 

F wave = P<V 2 T }V pTop , (3) 

where p is the density, F prop is the propagation velocity, and 
<Fj> is the root mean square velocity «F|> = 3/2^c)- If 
the wave is a sound wave, V prop is simply the sound speed, 
which is about 46 km s - 1 at the line formation temperature 
of C iv. For a plasma with solar abundances, hydrogen 


accounts for about 70% of the total mass, and at transition 
region temperatures n H & 0.8 n e (McWhirter, Thoneman, & 
Wilson 1975). Thus, the density p can be related to the 
electron density, n e , by the expression p = 1.14 m p n e , where 
m p is the proton mass. Noting that the electron pressure at 
10 5 K (in units of cm -3 K) is P e = 10 5 n e , we can then 
rewrite equation (3) as 

F s — 19m p P e £n C • (4) 

If the wave is an Alfven wave, F prop is the Alfven velocity, 
and equation (3) becomes 

F a = 1.4 x 10“ 3 (m p P b ) 1/2 ^nc & • (5) 

In Kolmogorov turbulent dissipation, energy cascades 
from the largest turbulent eddies to progressively smaller 
and smaller eddies until the bulk kinetic energy finally is 
converted to thermal energy through viscosity or Joule 
heating. The heating rate for the conversion of turbulent 
energy to thermal energy is (Shu 1992; Narain & Ulmsch- 
neider 1990, 1996) 

c, = ^, (6) 

where L is the scale size of the largest turbulent eddies (i.e., 
the scale size in which the turbulent energy first appears). A 
surface flux for turbulent energy can be defined as F turb = 
€ t h, where h is the thickness of the emission layer. We could 
estimate a value for h from plane-parallel atmospheric 
models of the Sun and other stars ( h ~ 10 km for the solar 
transition region lines; Nicolas et al. 1979), but it is very 
hard to know what value we should assume for L. Thus, we 
simply assume L = h, which corresponds to the situation in 
which the largest turbulent eddies span the width of the 
entire emission layer. This is a crude estimate, but it may 
not be unreasonable for the very thin transition region. 
Replacing p with l.\4m p n e , <F^> with 3/2^ c , and n e with 
10 ~ 5 P e , we then have 

F turb = 2.1 x 10- 5 m p E e ^c- (?) 

We now wish to compare F$, F A , and F turb with F rad , the 
radiative flux of the narrow component. The calculations of 
McWhirter et al. (1975) show that 

F rad = 2.4 x 10~ 22 E m (8) 

at the line formation temperature of C i v, T — 10 5 K, where 
E m is the emission measure. (For a homogeneous plane- 
parallel atmosphere, E m — nl h.) The observed C iv surface 
flux F c iv, which is computed assuming the C iv flux is 
emitted uniformly over the stellar surface, is related to the 
emission measure by the equation (Brown et al. 1984) 

F CIV = 1.65 x 10“ 23 E m a , (9) 

where a is the surface filling factor (0 < a < 1). Combining 
relations (8) and (9), we then have 

aF rad = 14.5F CIV (1 -/ BC // t ot) , (10) 

where the factor (1 —f B c/ftot) ensures that the C iv surface 
flux considered is for the narrow component only. 

Waves propagating through the transition region may 
deposit much of their energy in the corona as well as in the 
transition region, so we also wish to compare F s and F A 
with the coronal radiative flux, F cor , computed from the F x 
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values in Table 8. Based on the plasma models of Mewe, 
Gronenschild, & van den Oord (1985) and Mewe, Lemen, & 
van den Oord (1986), 

(6 x 10~ 23 £ m for log T cor = 6.0 
cor (2 x 10 -23 E m for log T cor = 7.0 , 1 j 

where T cor is the coronal temperature. For X-ray fluxes mea- 
sured in the Einstein IPC or ROSAT PSPC bandpasses, like 
those listed in Table 8, 

F x = 1.4x 10- 23 £ M « cor ^ (12) 

for both log T cor « 6.0 and log T cot « 7.0, where a cor is the 
surface filling factor for coronal emission and rj is a geo- 
metrical factor between 0.5 and 1.0. If the coronal emission 
is emitted close to the stellar surface, half the radiation from 
the plasma is emitted into the star rather than into inter- 
stellar space, and rj — 0.5. If the corona is very extended, 
however, and the star intercepts very little of the coronal 
emission, then rj « 1.0. For the Sun, rj is much closer to 0.5 
than to 1.0. Although this may not be the case for the other 
stars, we simply assume r) = 0.5 and combine equations (11) 
and (12) to find that 


8.6F X for log T cor = 6.0 
2.9 F x for log T cor = 7.0 . 


(13) 


Electron densities have been estimated from GHRS 
spectra for seven of the stars in our sample using different 
diagnostics (Maran et al. 1994; Wood et al. 1996; Linsky et 
al. 1995; Wood 1996; Dere & Mason 1993). Since P e should 
be roughly constant in the transition region (see, e.g., 
Jordan et al. 1987), electron pressures computed from these 
densities should apply at T = 10 5 K. Note that the lines 
used as density diagnostics have not been separated into 
narrow and broad components (the O iv] lines, for example, 
are too weak and/or too blended; see, e.g., Linsky et al. 
1995), and it is hard to judge how accurate densities com- 
puted from total line fluxes and ratios are in the presence of 
substantial inhomogeneity (see Judge, Hubeny, & Brown 
1997). 

Nevertheless, the electron pressures are listed in Table 9, 
along with values for aF rad , a cor F cor , F s , FJB, and F turb 
computed using equations (10), (13), (4), (5), and (7), respec- 
tively. In computing these quantities, we use the values of 
f civ, F xJ'Bc/fov and £ NC given in Table 8. We use the C iv 
values listed for / BC // to t> with the exception of a Cen A and B 
for which we are forced to use the Si rv values. For Capella, 
we assume that £ NC is equal to the sound speed (46 km s -1 ) 
rather than the velocity listed in Table 8, which we believe is 
an overestimate. We have to assume a value for T cor to 


TABLE 9 


Comparing Radiative, Wave, and Turbulent Fluxes 


Star 

log P e 
(cm- 3 K) 


Fluxes (10 s ergs cm 2 s 1 ) 



a co , F COI 

Fs 

FJB 

Ft urb 

AU Mic ... 

16.03 

16 

1000 

33 

4.4 

14 

HR 1099... 

14.99 

34 

730 

30 

13 

39 

Capella 

15.04 

15 

7.0 

30 

13 

36 

Procyon ... 

<14.92 

6.5 

14 

<23 

<11 

<27 

a Cen B . . . 

14.64 

0.67 

5.6 

3.0 

2.0 

1.8 

Sun 

15.00 

1.4 

2.9 

8.9 

3.9 

6.2 

a Cen A . . . 

14.79 

0.58 

1.0 

6.2 

3.4 

4.6 


estimate a cor F cor . For AU Mic, HR 1099, and Capella, we 
assume log T cor = 7.0 (Maran et al. 1994; Dempsey et al. 
1993b), and for the less active stars, we assume log T cor = 
6.0 (Mewe et al. 1995; Schrijver et al. 1995). 

Based on Table 9, for AU Mic and HR 1099 acoustic 
waves might barely be able to account for the narrow com- 
ponent radiative fluxes observed from the C iv emission 
layer, but they clearly cannot be responsible for even a 
significant fraction of the observed coronal emission. Thus, 
for AU Mic and HR 1099, acoustic waves are not a very 
promising source of heating. The situation is more ambigu- 
ous for the five less active stars listed in Table 9. For these 
stars, the estimated acoustic wave fluxes are large enough to 
account for much of the transition region and coronal radi- 
ation. Jordan et al. (1987) obtained similar results in their 
analysis of five main-sequence stars, including a Cen A and 
B. However, they pointed out that while acoustic waves 
may conceivably be able to account for the radiative losses 
from the stellar coronae, they cannot also account for 
expected conductive losses. More detailed analyses of the 
solar chromosphere and transition region using spatial- and 
time-resolved spectra have also cast doubt on the viability 
of the acoustic wave heating hypothesis (see, e.g., Athay & 
White 1978; Bruner 1978). We note that the F s values in 
Table 9 may actually be significantly lower than the F rad 
and F cor values for all the stars listed in the table if the 
surface filling factors, a and a cor , are significantly less than 1, 
which the solar example suggests is probably the case for 
the less active stars (especially for a C0T ). 

Another argument against the acoustic wave heating 
hypothesis is that we believe the magnetic field must play 
some role in increasing the narrow-component fluxes for 
more active stars. Although the narrow-component contri- 
bution to the C iv and Si iv lines decreases with increasing 
activity, the surface flux of the narrow component is still 
much larger for the active stars than for the inactive stars. 
For example, the narrow-component C iv surface flux of ft 
Cet is 6 times larger than the total C iv surface flux of ft 
Gem. It is difficult to explain this behavior without some 
type of magnetic heating process. Similar arguments against 
acoustic wave heating have been made in the past. The 
discovery that coronal activity is strongly correlated with 
stellar rotation clearly points toward a magnetic heating 
mechanism (see, e.g., Pallavicini et al. 1981). 

Alfven waves are a more plausible candidate for the 
narrow-component heating mechanism than acoustic 
waves. Table 9 demonstrates that it takes a magnetic field of 
only a few gauss for F A to exceed F s , and we can increase B 
even further in order to make F A comparable to F rad and 
F cor . For the Sun, a magnetic field of B « 10 G may be 
appropriate (Bruner 1978). For very active stars like AU 
Mic, photospheric magnetic fields as high as B & 4000 G 
have been found to cover a large fraction of the stellar 
surface (Johns-Krull & Valenti 1996). Jordan et al. (1987) 
also concluded that Alfven waves are a viable source of 
transition region and coronal heating, based on calcu- 
lations similar to ours. 

There are many other types of MHD waves besides 
Alfven waves that have been proposed as potential heating 
sources (see reviews by Narain & Ulmschneider 1990, 1996). 
In stellar atmospheres in which the Alfven speed can be 
expected to exceed the sound speed, slow-mode MHD 
waves and longitudinal wave modes in flux tubes do not 
propagate faster than the sound speed and therefore cannot 
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carry more flux than acoustic waves. Also, increasing the 
magnetic field does not increase the flux. Thus, all of the 
arguments presented above against acoustic wave heating 
also apply to heating by slow-mode and longitudinal MHD 
waves. In contrast, fast-mode MHD waves propagate at or 
above the Alfven speed, meaning that these waves, like the 
Alfven waves, are viable sources of transition region and 
coronal heating. 

A potential problem with wave-heating mechanisms is 
that one must be able to explain how the wave energy is 
dissipated. The usual picture is that after the waves are 
generated by subphotospheric convection, they travel 
upward through the atmosphere and eventually steepen 
into shock waves, thereby converting their energy to 
thermal energy through shock dissipation. However, the 
details of how this occurs have not yet been completely 
worked out. 

In contrast, if the narrow component nonthermal veloci- 
ties are interpreted as turbulent velocities, conversion of the 
turbulent energy to thermal energy is not difficult to 
explain. One would expect turbulence to be naturally dissi- 
pated through a Kolmogorov cascade process. The values 
of F turb in Table 9 suggest turbulent heating may indeed be 
a viable heating mechanism, since as long as a is not too 
small, F turb « F rad . There is no point in comparing F turb 
with F cor , because turbulent energy should be dissipated 
locally. Thus, for the turbulent dissipation interpretation to 
work, we have to show only that F turb ~ F rad , in contrast to 
the wave fluxes, for which we want F wave > F rad + F cor . 
Although F turb does appear to be comparable to F fad , we 
caution that uncertainties in the calculation are high, as 
equation (7) is more of an order-of-magnitude kind of equa- 
tion than equations (4) and (5). In a similar calculation, 
Nicolas et al. (1979) also found that turbulent dissipation 
may be able to account for the transition region heating on 
the Sun, assuming once again that a is not too small. We 
now present another argument in favor of the turbulent 
heating hypothesis. 

Jordan et al. (1987) found empirically that 

( 14 ) 

for main-sequence stars. Combining this relation with the 
relation 

fcvOcFi,;,]* 012 (15) 

from Ayres et al. (1995), we find that 

F ciW ozPl- 06±0 - 01 g- x - 11±0A2 . (16) 

Since equation (2) provides us with an empirical relation 
between £ NC and g, we can express equation (7) in terms of 
P e and g. If we then equate F turb and F rad in equations (7) 
and (10), we can derive an expression analogous to equation 
(16) and test the F rad = F turb assumption. Since equation 
(16) originates from analyses that assume homogeneity and 
spherical symmetry, for the purposes of this comparison we 
assume a = 1 and (1 — / BC //tot) = 1. We then have 

F cl yOcP e g^ ±0 - 21 . (17) 

Equations (16) and (17) agree very well, which provides 

support for the F rad — F turb assumption. It is possible that 
the narrow-component heating mechanism, whatever it 
may be, deposits its energy into the transition region 
through turbulence, and the transition region is then heated 
by the dissipation of this turbulence. We are not sure what 


the source of the turbulent energy is, but as mentioned 
above, we believe that magnetic activity must be involved. 
Considering the vast literature that exists concerning pos- 
sible coronal heating mechanisms, it should come as no 
surprise that there are many theories that are compatible 
with the requirements that the narrow-component heating 
mechanism be connected with magnetic activity and that 
the energy first be deposited into the upper atmosphere in 
the form of turbulence. For example, Hollweg (1984) pro- 
posed that Alfven waves are dissipated in the corona 
through Kolmogorov turbulence. Also, Gomez & Ferro 
Fontan (1988) proposed that MHD turbulence can be pro- 
duced by motions of magnetic loop footpoints induced by 
subphotospheric convection. 

The empirical £ N<? ocg~ 0 ' 68 relation that we have 
observed may be an important diagnostic for the narrow 
component heating mechanism. In the turbulent dissipation 
scenario, it is unclear why this relation should exist. Any 
proposal for the source of the turbulent energy should seek 
to explain the apparent gravity dependence. In the wave- 
heating scenario, the gravity dependence can be understood 
through the probable source of the waves — 

subphotospheric convection, because the convective energy 
should be dependent on the surface gravity (see, e.g., Stein 
1981 ; Ulmschneider & Stein 1982). 

6. SUMMARY 

Although the main intent of this paper is to provide a 
global comparison of the transition region line profiles of 
late-type stars, for a few stars we do present new results 
specific to those stars: 

1. We find that the emission-line redshifts of a Cen A and 
B increase with line formation temperature up to 
log T = 5.2. These redshifts are almost identical to those 
previously observed for the Sun and Procyon, although the 
redshifts of a Cen B are about 4 km s -1 larger than 
the redshifts of the other three stars. Above log T = 5.2, the 
redshifts decrease dramatically for the Sun, a Cen A, and a 
Cen B, but for Procyon, they continue to increase. 

2. Using the O rv] lines, we derive electron densities at 
log T = 5.14 of log n e = 9.65 ± 0.20 and log n e = 9.50 
± 0.30 for a Cen A and a Cen B, respectively, correspond- 
ing to pressures (in units of cm -3 K) of log P e — 14.79 
± 0.20 and log P e = 14.64 ± 0.30. These pressures are con- 
sistent with coronal pressures measured from EUVE data, 
contrary to what has been found for Capella and Procyon. 

3. GHRS observations of the C iv lines of ^ Dra 
obtained in three different years (1992, 1993, and 1995) 
show substantial line profile variability. The reason for this 
variability in uncertain, but rotational modulation is a 
possibility. 

By comparing the Si iv and C rv profiles of many stars, 
we have conclusively shown that broad line wings are mani- 
festations of stellar activity, as the prominence of these 
wings is clearly correlated with both the C iv and X-ray 
surface fluxes. We have tried to quantify the generally non- 
Gaussian Si iv and C iv profiles with two-Gaussian fits. We 
believe that microflaring is the best interpretation for the 
broad components of these two-Gaussian fits, for the fol- 
lowing reasons : 

1. The very large widths of the broad components 
suggest supersonic velocities. We are reluctant to interpret 



764 


WOOD, LINSKY, & AYRES 


Vol. 478 


these widths as being due to supersonic turbulence, because 
supersonic turbulence cannot exist for long due to shock 
dissipation. By associating the broad components with 
microflaring, we can explain the broad lines as being due to 
jets of plasma expelled at high velocities from magnetic 
reconnection sites. 

2. If microflares are in fact common in stellar transition 
regions, it would not be surprising if they were especially 
prevalent in the transition regions of very active stars. Thus, 
the correlation we observe between the broad components 
and stellar activity is consistent with this interpretation. 

3. Solar transition region explosive events, which can be 
considered a form of microflaring, might be the solar analog 
for the broad components. For the solar-like star a Cen A, 
the broad component has a width of 109 ± 10 km s _1 and 
is blueshifted relative to the narrow component by 9 + 3 
km s^ 1 . Both of these properties are consistent with the 
average properties of the solar explosive events. However, 
for a Cen A, the broad component accounts for about 
25% + 4% of the total line emission, while it is currently 
thought that the explosive events account for only about 
5% of the solar Si iv and C iv emission. Perhaps two- 
Gaussian fits overestimate the broad-component contribu- 
tion to the line flux and/or perhaps the explosive event 
contribution to the Sun’s transition region emission has 
been underestimated due to insufficient spatial and time 
resolution. 

4. In addition to the solar explosive events and the broad 
components, other evidence supports the presence of con- 
tinuous microflaring in the atmospheres of cool stars. For 
example, very active stars typically have coronal tem- 
peratures above 10 7 K, but on the Sun temperatures this 
high are seen only during flares (see, e.g., Haisch, Strong, & 
Rodono 1991). Observed correlations between microwave 
and X-ray emission are similar for both active stars and 
solar flares, which is also suggestive of microflaring (Giidel 
1994; Benz & Giidel 1994). 

Finally, we analyze the properties of the narrow com- 
ponents of the two-Gaussian fits and discuss possible 
heating mechanisms for the narrow component emission. 
Our main findings are as follows : 
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1. When we convert the widths of the narrow com- 
ponents to nonthermai velocities (£ NC ), we find that the 
resulting velocities are correlated with stellar surface gravity 
(i g ). For the low-gravity evolved stars, our data are consis- 
tent with constant £ NC equal to the sound speed («46 km 
s“ *). For the high-gravity main-sequence stars, £ NC is less 
than the sound speed and £ NC oc g~ 0 68 . 

2. Using the empirical relation between £ NC and g, we can 
reproduce scaling laws computed for main-sequence stars 
by Jordan et al. (1987) from atmospheric models based on 
IUE data, if we assume that the radiative flux is equal to the 
turbulent flux. This suggests that the narrow components 
may be associated with heating through turbulent dissi- 
pation. Order-of-magnitude calculations suggest that the 
turbulent and radiative fluxes are indeed roughly equal for 
the seven stars in our sample that have measured transition 
region densities, as long as the surface filling factor of the 
transition region emission is not too low. 

3. If the nonthermai velocities, £ NC , are due to motions 
induced by Alfven waves rather than turbulence, it is quite 
possible that the waves carry enough flux to account for the 
narrow-component emission and much of the coronal emis- 
sion. Thus, Alfven wave dissipation is another viable 
heating mechanism. 

4. Results in this paper and in previous work suggest that 
acoustic waves are not a likely source of the narrow- 
component emission. For one thing, our data demonstrate 
that stars with different activity levels that are otherwise 
very similar can have significantly different narrow com- 
ponent fluxes, which suggests that the heating mechanism is 
probably associated with magnetic activity. 

An important goal for future work is to pursue these issues, 
particularly for the Sun now that SOHO is available for 
high-resolution spectroscopic studies of our nearby star. 
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